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ABSTRACT 


A  Horizontal  Line  Array  (HLA)  receiver  is  subject  to  left-right  bearing  ambiguity.  It  is 
also  susceptible  to  bearing  bias  caused  by  non-zero  vertical  arrival  angles.  The  true  bearing 
angle  and  the  vertical  arrival  angle  may  both  be  estimated  if  a  second  observation  is  made, 
with  the  HLA  rotated  to  a  new  orientation.  A  closed-form  solution  is  presented  for  the  case 
of  stationary  source  and  receiver.  For  small  observation  errors,  the  resulting  errors  in  the 
estimated  angles  can  also  be  expressed  in  closed  form. 


r£sum£ 

Les  recepteurs  de  sonar  en  ligne  horizontale  (Horizontal  Line  Array,  HLA)  sont  sujets  a  des 
ambiguit4s  a  gauche  ou  4  droite  quant  au  relevement.  Ds  sont  £galement  sujets  a  des  erreurs  de 
relevement  causees  par  des  angles  d’arrivee  non  nuls  par  rapport  a  la  verticale.  Le  relevement 
veritable  et  Tangle  d’arrivee  par  rapport  a  la  verticale  peuvent  tous  deux  etre  estimes  si  on 
effectue  une  deuxieme  observation  apres  avoir  aligne  le  HLA  Bur  un  nouveau  cap.  Une  solution 
analytique  complete  est  pr6sent4e  pour  le  cas  d’une  source  et  d’un  recepteur  stationnaires.  Pour 
de  petites  erreurs  d ’observation,  les  erreurs  r£sultantes  des  angles  estimes  peuvent  6galement 
etre  exprimees  sous  forme  de  solution  analytique  complete. 
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1  Introduction 


A  Horizontal  Line  Array  (HLA)  receiver  is  a  one-dimensional  array  of  sensors  lying  in  the 
horizontal  plane.  An  HLA  can  be  used  to  detect  an  underwater  acoustic  source,  and  to  obtain 
an  estimate  of  the  direction  of  the  arriving  signal.  Due  to  the  one-dimensional  nature  of  the 
HLA,  the  bearing  estimates  are  subject  to  left-right  ambiguity.  HLA  bearing  estimates  are 
also  susceptible  to  bearing  bias  caused  by  non-zero  vertical  arrival  angles. 

The  true  bearing  angle  and  the  vertical  arrival  angle  may  both  be  estimated  if  a  second 
observation  is  made.  This  can  be  accomplished  by  rotating  the  HLA  to  a  new  orientation.  The 
problem  can  be  modelled  as  two  nonlinear  equations  in  two  unknowns.  A  closed-form  solution 
is  presented  for  the  case  of  stationary  source  and  receiver.  Error  bounds  on  the  estimates 
can  be  calculated  as  a  function  of  observation  errors,  HLA  rotation  angle,  and  signal  arrival 
geometry.  For  small  observation  errors,  the  resulting  errors  in  the  estimated  angles  can  aiso 
be  expressed  in  closed  form.  — 


2  The  HLA  Bearing  Bias  Problem 


Figure  2.1  illustrates  the  geometry  of  a  plane  wave  signal  arriving  at  a  Horizontal  Line 
Array.  The  HLA  is  lying  in  the  plane  *  =  0.  The  azimuthal  angle  0  is  measured  in  the  plane, 
clockwise  from  the  HLA  axis.  The  arriving  signal  has  azimuthal  angle  9  and  vertical  arrival 
angle  <f>. 

The  HLA  sensors  may  be  processed  using  one  of  several  methods,  to  provide  an  esti¬ 
mate  of  cos/9,  where  0  is  the  angle  of  incidence  between  the  array  axis  and  the  signal  vector 
[Walker,  1985].  The  angle  of  incidence  is  a  function  of  the  azimuthal  angle  and  the  vertical 
angle: 


cos  0  =  cos  6  cos  <f>. 

(2.1) 

The  estimate  of  arrival  direction  therefore  has  a  left-right  ambiguity: 

0  =  ±0o 

(2.2) 

where 

0o  =  cos-1  (costfcos*)  ;  0  <  0o  <  180°  . 

(2.3) 

If  the  signal  vector  is  in  the  plane  <j>  =  0,  then  cos/9  =  cos  9.  However,  if  the  vertical 
arrival  angle  <f>  is  non-zero,  then  the  measured  angle  of  incidence  0  may  not  be  equal  to  the 
azimuthal  angle  9.  0  is  therefore  a  biased  estimator  of  9.  The  effect  of  non-zero  <f>  is  to  make 
cos  0  smaller  than  cos  0.  Therefore  0  is  biased  toward  the  broadside  direction  cos  0  =  0. 

A  graph  will  be  used  to  show  the  magnitude  of  the  bias,  for  9  between  0°and  180®.  The 
bias,  B,  is  defined  as  follows: 

B=  0-8 

=  cos_1(cos0cos*)  -  9  ;  0°  <  0  <  180°  (2.4) 

=  -  cos-l(cos0cos0)  -  0  ;  -180°  <  0  <  0°  . 

Figure  2.2  shows  the  bias  error  B  as  a  function  of  9,  for  several  values  of  the  vertical  arrival 
angle  <f>.  It  is  evident  from  this  figure  that  bearing  bias  can  be  a  serious  problem.  The 

maximum  bias  occurs  in  the  endfire  directions,  0  =  0°and  0  =  180°;  the  maximum  absolute 

bias  error  is  equal  to  the  vertical  arrival  angle  |*j: 

I B\  <  |*|  •  (2.5) 
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BIAS  B  IB  DEGREES 


Figure  2.1:  Arrival  Angle  Geometry 


AZIMUTHAL  ANGLE  0  IN  DEGREES 


Figure  2.2:  HLA  Bearing  Bias  as  a  function  of  azimuthal  angle  0,  for  several  values  of  the 
vertical  arrival  angle  6- 


3  Estimation  of  Azimuthal  and  Vertical 
Arrival  Angles 


3.1  Solution  by  Rotation  of  the  Array 

The  true  azimuthal  angle  and  the  vertical  arrival  angle  may  both  be  estimated  if  a  second 
observation  is  made.  This  can  be  accomplished  by  rotating  the  Horizontal  Line  Array  to 
a  new  orientation.  The  problem  may  then  be  approximated  as  two  nonlinear  equations  in 
two  unknowns.  A  closed-form  solution  can  be  obtained  for  the  case  of  stationary  source  and 
receiver.  For  small  errors  in  the  observations  and  in  the  assumed  rotation  angle,  the  resulting 
errors  in  the  estimated  angles  can  also  be  expressed  in  closed  form. 


3.2  Assumptions  and  Restrictions 

A  number  of  simplifying  assumptions  are  made  in  order  to  find  a  computable  closed-form 
solution: 

(a)  The  array  is  assumed  to  be  perfectly  stable;  errors  due  to  tilting  and  stretching  of  the 
array  are  not  included.  However,  the  standard  deviation  of  the  angular  observation  error 
may  be  increased  to  approximate  this  effect. 

(b)  The  source  and  receiver  positions  are  assumed  to  be  stationary.  However,  the  effect  of 
minor  relative  motion  may  be  approximated  by  enlarging  the  error  bounds  for  the  final 
estimated  angles. 

(c)  The  propagation  path  is  assumed  to  be  the  same  for  both  observations;  that  is,  the  vertical 
arrival  angle  <f>  remains  constant.  This  assumption  restricts  the  usefulness  of  the  method 
to  environments  in  which  the  vertical  arrival  angle  is  stable. 

(d)  The  rotation  angle  AO  is  assumed  to  be  larger  than  the  magnitude  of  the  observation 
errors  e,  i.e.  |sin  A0|  »  e. 
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3.3  The  Analytic  Model 

The  Horizontal  Line  Array  is  rotated  through  an  angle  A 9,  where  sinA0  ^  0.  The  ob¬ 
servations  cos  fa  and  cos  fa  are  made  before  and  after  rotation,  respectively.  Let  9\  and  6 j  be 
the  azimuthal  angles  clockwise  from  the  array  forward  endfire  direction  before  and  after  rota¬ 
tion.  Let  4> l  and  fa  be  the  corresponding  vertical  arrival  angles  before  and  after  rotation.  Then 


it  is  possible  to  describe  the  problem  by  two  equations  in  four  unknowns: 

cos/Sj  =  cosfljcos^i  (3.1) 

cos  fa  =  cos  9-i  cos  <f> 2  .  (3.2) 

The  number  of  unknowns  may  be  reduced  by  assuming  fa  =  fa  =  <f>  and  9 1  =  9,  9i  = 
9  -  A0,  to  obtain  two  equations  in  two  unknowns: 

cos  fa  =  cos  9  cos  <f>  (3.3) 

cos  fa  —  cos  (9  -  A 9)  cos  <j>.  (3.4) 


The  assumption  of  constant  vertical  arrival  angle  <f>  is  reasonable  if  the  measurements 
are  made  simultaneously  by  a  pair  of  collocated  HLAs,  or  if  the  acoustic  propagation  path  is 
stable.  For  instance,  in  underwater  acoustics,  bottom  reflected  paths  arriving  at  steep  vertical 
angles  are  relatively  stable. 

The  assumpt'on  of  constant  azimuthal  angle  9  is  reasonable  if  there  is  little  relative  move¬ 
ment  of  source  and  receiver  between  observations.  That  is, 

(a)  the  measurements  are  made  simultaneously  by  a  pair  of  collocated  HLAs,  or 

(b)  the  source  and  receiver  are  stationary,  or 

(c)  the  source  is  in  the  far  field. 

In  practice,  a  finite  amount  of  time  At  would  be  required  to  rotate  the  array.  The  array  may 
also  be  shifted  a  distance  D.  Therefore  an  error  Ee  may  be  introduced  in  the  estimate  of  6 
due  to  motion  of  a  source  with  velocity  V  at  range  R: 

_  <  180  fD  +  VAt\  . 

£>  “  T  (— R— )  P  5> 

This  error  should  be  included  when  calculating  the  total  error  in  the  final  estimate  of  9.  For 
example,  if  D  =  0.075nm,  V  =  llkt,  At  =  O.lhr  and  R  =  3nm,  then  E$  «  3.3°. 


$ 
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3.4  The  Analytic  Solution 

The  analytic  solution  may  be  found  by  expanding  equation  3.4  to  obtain: 


cos  0  cos  4>  =  cos  f}\ 


sin  0  cos  4>  = 


cos  Pi  -  cos  Pi  cos  A0 
sin  A  6 


The  vertical  arrival  angle  <f>  is  found  by  squaring  and  summing  equations  3.6  and  3.7  to  obtain: 


,  /  „  /  cos  Pi  —  cos  p\  cos  A 0 

m *-  (“tfM  - 


If  cos  <f>  ^  0,  then  the  azimuthal  angle  0  is  uniquely  defined  by: 


!  0  <  <f>  <  90 


COS  0  —  COS  Pi/  cos^  -39. 

sin0  =  (cos  pi  —  cos  P\  cos  A0)/ (sin  A0  cos  '  '  ' 

If  <f>  =  90°,  then  0  is  indeterminate.  Also,  when  observation  errors  are  present,  care  must  be 
taken  to  ensure  that  |cos^|,  |  cos  8\  and  |sin0|  do  not  exceed  unity. 

The  next  section  shows  graphs  of  the  relationship  between  Pi,  Pi  and  <f>,  and  illustrates 
the  sensitivity  of  the  solution  to  problem  geometry. 

Section  5  returns  to  the  analytic  solution  to  derive  closed-form  expressions  for  maximum 
and  r.m.s.  errors  in  the  estimated  values  of  azimuthal  angle  0  and  the  cosine  of  the  vertical 
angle  4>. 
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4  Graphic  Representation 


4.1  Graphic  Estimation  of  Vertical  Arrival  Angle 

The  HLA  angle  estimation  problem  can  be  solved  directly  using  a  computer.  However,  it 
is  instructive  to  examine  graphs  of  the  relationships  between  the  parameters  of  interest.  Such 
graphs  provide  an  appreciation  for  the  sensitivity  of  the  solution  to  observation  errors.  Graphs 
may  also  provide  an  alternative  to  using  a  computer.  However,  the  use  of  such  graphs  may  be 
cumbersome  due  to  the  large  number  of  variables  [fa,  fa,  AO ,  6,  and  <f>). 

To  obtain  a  graphic  solution,  it  is  convenient  to  first  estimate  the  vertical  angle  4>.  There 
are  several  ways  to  solve  graphically  for  <f>.  Four  methods  will  be  considered  here.  In  each 
of  these  methods,  separate  graphs  Me  required  corresponding  to  discrete  values  of  the  rotation 
angle  AO. 

Method  (a)  In  this  method,  is  plotted  against  fa  for  several  discrete  values  of  <f>.  This 
allows  graphic  resolution  of  left-right  ambiguity,  but  requires  a  quadrant  search  for  the 
correct  solution. 

Method  (b)  This  method  uses  absolute  angles  only.  |£i|  is  plotted  against  \fa\  for  several 
discrete  values  of  <f>.  This  graph  may  also  be  used  to  resolve  left-right  amibiguity,  but 
requires  careful  interpretation.  The  larger  scale  is  easier  to  read. 

Method  (c)  This  method  uses  a  nonlinear  scale,  cos  fa  is  plotted  against  cos  fa  for  several 
discrete  values  of  <f>.  This  is  similar  to  Method  (b),  except  for  the  use  of  the  cosine 
scale.  Horizontal  Line  Array  resolution  capability  is  usually  uniform  in  the  cosine  do¬ 
main.  That  is,  cos  f)  can  be  estimated  with  an  accuracy  which  is  independent  of  0.  The 
graph  for  Method  (c)  therefore  provides  a  better  feeling  for  sensitivity  to  observation 
errors. 

Method  (d)  This  method  is  based  on  the  deviation  between  the  expected  and  observed  values 
of  |/3j|.  Given  the  observation  |/?i|,  two  forecasts  are  made  for  |/?j|,  for  right  and  left 
bearing  hypotheses,  under  the  assumption  that  $  =  0.  After  rotation,  the  forecast 
which  agrees  most  closely  with  the  observed  value  \fa\  determines  the  choice  of  right  or 
left  bearing.  If  the  agreement  is  not  exact,  the  difference  may  be  used  to  estimate  the 
magnitude  of  <f>.  This  method  allows  large  values  of  <j>  to  be  identified. 


4.1.1  Example  of  Method  (a) 


Figure  4.1  shows  the  graphs  for  Method  (a),  for  rotation  angles  of  20°,  45°,  and  90°.  The 
observed  values  of  cos  0i  and  cos  02  do  not  immediately  indicate  the  signs  of  these  angles.  How¬ 
ever,  for  any  pair  of  observations,  there  is  only  one  signed  pair  that  fits  the  data.  The  answer 
is  found  by  searching  through  the  four  possible  sign  pairs: 


(+iai,  +m 

(+!0il,  -|0z|) 

HAI,  +IAI) 

hai,  -m 

Consider  the  following  example: 

A9  =  20°,  |0X|  =  96.4°,  |02 1  =  108.8° 


In  figure  4.1,  the  only  solution  pair  is: 

0i  =  -96.4°,  02  =  -108.8° 

which  lies  on  the  contour  <f>  —  50°.  These  values  are  difficult  to  read  unless  the  plot  is  quite 
large.  The  estimated  azimuthal  angles,  from  Figure  4.5,  are: 

Bi  =  -100°,  =  -120°. 

4.1.2  Example  of  Method  (b) 

Figure  4.2  shows  the  corresponding  graphs  for  Method  (b).  Given  |0x|  and  02|  as  before, 
the  vertical  angle  <f>  may  be  read  from  the  graph  by  interpolating  between  contours.  This  graph 
is  easier  to  read  because  the  contours  are  continuous  and  are  plotted  on  a  larger  scale.  The 
left-right  ambiguity  is  resolved  as  follows.  For  any  given  vertical  angle  <f>,  there  are  two  possible 
a  priori  values  of  |02l  f°r  each  observed  angle  0\.  The  largest  value  corresponds  to  rotation 
away  from  the  signal  direction  {9\  and  A9  have  opposite  signs);  the  smallest  value  corresponds 
to  rotation  towards  the  signal  direction  (#i  and  A 9  have  the  same  sign). 


4.1.3  Example  of  Method  (c) 

Figure  4.3  shows  the  graphs  for  Method  (c).  The  vertical  angle  4>  and  the  left-right 
ambiguity  may  be  solved  as  for  Method  (b).  The  natural  cosine  scale  gives  a  better  appreciation 
of  the  ability  to  resolve  different  values  of  For  instance,  a  large  value  of  <f>  produces  a  large 
shift  in  the  observed  value  0\  when  the  signal  is  arriving  from  the  endfire  directions  9  =  0°and 
9  =  180°;  however  the  HLA  is  less  capable  of  resolving  these  differences  in  the  endfire  directions. 
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Figure  4.3:  Estimation  of  vertical  arrival  angle  <j>  using  graphic  Method  (c) 


4.1.4  Example  of  Method  (d) 

Figure  4.4  shows  the  graphs  for  Method  (d).  In  this  Method,  the  a  priori  estimate  |02|,  is 
based  on  the  assumption  that  the  vertical  arrival  angle  ^  is  equal  to  0.  There  are  two  possible 
solutions,  corresponding  to  the  ambiguous  left  and  right  solutions: 

|0z|(a)  =  cos'1  (cos(|£jj  +  |A0|))  ;  =>  faAO  <  0  (4.1) 

or 

IA|(t)  =  cos-1  (cos(|0j|  —  |A0|)))  ;  => &A0  >  0  .  (4.2) 

After  the  observation  |02|,  the  chosen  estimate  |02|  is  the  one  which  agrees  most  closely  with 
the  observation.  The  deviation  6  is  then  defined  as  follows: 

*  =  m  -  \h\  (4.3) 

Note,  in  Figure  4.2,  that  zero  deviation  is  expected  when  <f>  =  0.  If  4>  ^  0,  then  a  negative 
deviation  is  expected  when  the  rotation  is  away  from  the  signal  direction;  a  positive  deviation 
is  expected  when  the  rotation  is  towards  the  signal  direction. 

Consider  the  same  example  again: 


AO  =  20°,  |0i|  =  96.4°  and  |02|  =  108.8° 


Then: 


|&|(o)  =  116.4°  and  |/32|(6)  =  76.4° 


The  appropriate  choice  is 


\Pi\  =  \fo\[a)  =  116.4°, 


which  implies  that  0\AO  <  0  and  therefore 

01  =  -96.4°. 

However,  there  is  not  a  perfect  match  between  the  predicted  angle  before  rotation  and  the 
observed  angle  after  rotation.  The  deviation  is: 

6  =  |108.8°j  -  |116.4°|  =  -7.6°. 

From  Figure  4.4,  the  value  of  the  vertical  arrival  angle  4>  is  therefore  50°.  Note  that  in  Fig¬ 
ure  4.4,  the  derivative  discontinuties  correspond  to  the  corners  of  the  rectangles  in  Figure  4.2. 
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4.2  Graphic  Estimation  of  Azimuthal  Angle 


Once  the  value  of  the  vertical  arrival  angle  ^  has  been  estimated,  it  may  then  be  used  to 
correct  0  to  estimate  the  true  bearing  9.  Equation  2.1  defines  the  relationship  between  6,  4> 
and  0.  Figure  4.5  shows  the  relationship  between  9  and  0  for  several  values  of  4>. 

To  continue  the  example  from  the  preceding  section,  it  has  been  shown  that,  if 
A 6  =  20°,  |  =  96.4°,  and  \02\  =  108.8°, 

then  it  can  be  deduced  that 

0i  =  -96.4°,  02  =  -108.8°,  and  <f>  =  50°  . 

Equation  2.1  or  Figure  4.5  may  be  used  to  obtain  the  estimated  azimuthal  angles: 

9\  =  -100°,  e2  =  -120°  . 
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Figure  4.5:  The  true  bearing  9  can  be  obtained  from  observation  0,  if  the  vertical  angle  <t>  is 
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5  Error  Estimates 


5.1  Model  of  Observation  Errors 


The  analytic  solution  provides  an  exact  closed  form  for  the  azimuthal  and  vertical  arrival 
angles.  However  there  may  be  errors  in  the  observed  values  cos  0\  and  cos  02  as  well  as  error 
in  the  assumed  rotation  angle  A 6. 


5.1.1  Errors  in  cos  0 


It  will  be  assumed  that  the  standard  deviation  of  the  HLA  observation  cos/9  is  independent 


Q  co§  0  —  £  —  sin  ao 


where  a0  is  the  standard  deviation  of  the  bearing  error  in  the  broadside  direction  |/9|  =  90°,  for 
a  specified  signal-to-noise  ratio.  Figure  5.1  shows  the  observation  errors  A/9  corresponding  to 
cos/9  ±  t  as  a  function  jf  0  for  the  case  ao  =  1°-  Note  that  maximum  errors  |A^|max  occur  in 
the  endfire  directions,  and  also  at  angles  0+  |A/9|max  and  180  -  | A0|mmx.  For  small  values  of  e, 
a  second-order  Taylor  series  expansion  of  cos/9  leads  to  an  approximate  expression  for  |A/9jm»x: 


Acos£  es  —  sin  0  (  —  A/9 )  -  -  cos/9  [  —  A/9 ) 
\  180  /  2  \  180  / 


*  -  sin/9  +  signfsin  0)\J  { sin2  0  ±  2  cos  0 A  cos  0)  + 

—  A0  ss  - - -  •  cos/9  ^  0  (5.3) 

iso  cos  0  v  ' 


~|A£|max  «  >/2e 

ISO 


max  J — 2a0  c»  lOJy/ffo  • 


In  this  example,  the  maximum  error  is  approximately  10.7° 
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Figure  5.1:  Bounds  on  error  in  the  observed  angle  p  ,  corresponding  to  cos/9  ±  sinero,  where 

a0  -  r. 


5.1.2  Errors  in  Rotation  Angle  A 6 


The  standard  deviation  in  the  rotation  angle  A0  will  be  denoted  by 


ct&s  =  h 


where  h  is  given  in  degrees. 


5.2  Resulting  Errors  in  the  Closed-Form  Solution 


For  suitably  small  errors  t  and  h,  the  resulting  errors  in  the  estimated  azimuthal  and 
vertical  arrival  angles  0  and  ^  can  be  approximated  in  closed  form.  However,  if  the  errors  are 
large,  more  elaborate  methods  are  required  to  estimate  the  errors  [Theriault,  1987], 

First  derivatives  are  obtained  by  differentiating  equations  3.6  and  3.8  with  respect  to 
cos^i,  cos  fa  and  A 0  to  obtain  the  following  expressions: 


d  cos  <f>  sin  82 

d  cos  sin  A  8 

(5.7) 

d  cos  <f>  sin  8\ 

d  cos  P2  sin  A  8 

(5.8) 

d  cos  <f>  r  sin  8\  sin  82  cos  <p 

8A8  iso  sinA0 

(5.9) 

d$i  180  COS  82 

d  cos  Pi  r  cos  4>  sin  A 8 

(5.10) 

88i  180  COB 

d  cos 02  *  cos  ^  sin  A0 

(5.11) 

30j  cos  0i  sin  0 2 

d  A0  sin  A 0 

(5.12) 

The  error  analysis  is  based  on  a  first-order  Taylor  series  expansion  in  the  errors  t  and 
h.  It  is  assumed  that  these  errors  are  small  with  respect  to  the  rotation  angle: 


e  =  sin  (To  4C 


[  sin  A0|cosd> 
jsinA0|  . 


(5.13) 


It  is  further  assumed  that  the  errors  in  cos^i,  cos^i  and  the  error  in  the  assumed  rotation 
angle  A 0,  are  statistically  independent.  Two  cases  will  be  considered  here: —  maximum  errors 
and  root-mean-square  errors. 


5.2.1  Maximum  Errors 

It  will  be  assumed  that  e  and  h  are  strict  upper  bounds  on  the  errors  in  cos/?  and  A 0: 


(cos  P  -  cos^j  <  E 
|A0  -  A0|  <  h  . 


(5.14) 


Then  upper  bounds  on  the  absolute  errors  in  the  vertical  angle  cos  4>  and  the  azimuthal  angle 
9  are  given  by  the  first-order  Taylor  expansions: 


<  |  d  cos  4>  |  1  d  cos  4>  |  [  d  cos  4>  | 


I  A  ii  ''W  Y"  ,  v  tvo  u/ 

cos <f>  -  cos 4>\  cs  -= - -g-  e  +  -= - 3-  e  + 

1  0  cos  pi  d  cos  £2 


<  Pain#,!  -f  |sin  fl2l~|  sin  0,  sin  02  cos  <j> 

W  |sinA0|  C  180  sin  AS 


(5.15) 


-  aSi7'  + 


,  as, 

ai  acos^2  e+  aAS 


I  CQsSil  -f  I  COS  S2|  1  _T 
|  sin  AS{  cos  4>  181 


I  cos  S,  sin  S2  I 


180  sin  AS 


(5.16) 


=  tf#  • 

Annex  A  contains  a  parametric  analysis  of  maximum  errors. 

5.2.2  R.M.S.  Errors 

The  r.m.s.  errors  of  the  estimated  azimuthal  angle  and  the  cosine  of  the  vertical  angle 
will  be  denoted  by  S$  and  Scol$  respectively.  For  small  error  standard  deviations  e  and  h,  a 
first-order  Taylor  series  expansion  may  be  used  to  obtain  the  following  estimated  errors: 


Scot  6 


V  a  cos  a,  /  v  a  cos  a*  )  \  aAS  j 


(5.17) 


(sin2  S,  sin2  S2)e2  +  ( —  sin  S,  sin  S2  cos  4>)2h2 
„ _ 180 _ 

I  sin  ASI 


iso5*  w 


,a  cos  a. 


-),+(>s^),+(a*)T 


(cos*  0,  +  cos*  02)es  +  (^cos0,  sin  02  co8<f>)2h2 
|sin  A0|costf> 


(5.18) 


Note  that  the  approximate  error  in  cos  <£  is  not  exact  as  4>  approaches  zero  because  of 
the  boundary  constraint  |cos^j  <  1.  Similarly,  the  approximate  error  in  0  is  not  exact  as 


,  A*  .V  .  (A  A  AA 


.  a  «v -v a/..a/.  •-./ 


■  *  *  11  .  *  *  w 

ifovM 


0  approaches  the  endfire  directions.  In  these  situations,  the  approximate  error  expression 
overestimates  the  actual  r.m.s.  error. 


From  these  equations,  the  effect  of  an  error  h  in  the  rotation  angle  is  comparable  to  the 
effect  of  the  broadside  angular  resolution  <ro.  Therefore,  if  h  cos  <f>  >  <tq  then  the  error  in  rotation 
angle  dominates,  while  if  hcos#  <  ctq  the  angular  resolution  of  the  array  dominates.  Errors  in 
the  rotation  angle  are  less  critical  when  the  vertical  arrival  angle  is  steep. 

Annex  B  contains  a  parametric  analysis  of  r.m.s.  errors.  Annex  C  contains  examples  of 
the  effect  of  error  in  rotation  angle. 


5.3  An  Example  of  Maximum  Errors 

5.3.1  Errors  When  Vertical  Angle  Bias  is  Not  Removed 

First,  for  purposes  of  comparison,  it  is  instructive  to  consider  the  errors  which  occur  if  the 
array  is  not  rotated,  and  the  vertical  arrival  angle  <f>  is  assumed  to  be  zero.  For  the  special 
case  A 9  =  0,  the  maximum  errors  are  defined  as  follows: 


U,  *  \B(9,4)\  +  \A0{$,M\  +  h  ;  A9  =  0  (5.19) 

U#  =  <f>  ;  A9  =  0  (5.20) 

Figure  5.2  shows  the  approximate  error  bound  Ut  on  azimuthal  angle,  for  <f>  —  20°.  The 
broadside  resolution  angle  is  tr0  =  1°,  and  the  error  h  in  rotation  angle  is  assumed  to  be 
zero.  The  error  is  most  severe  in  the  endfire  directions.  At  endfire,  the  estimates  are  subject 
to  the  maximum  bias  error  plus  the  maximum  angular  resolution  and  rotation  angle  error. 

5.3.2  Errors  in  Estimated  Azimuthal  Angle  after  Rotation 

Figure  5.3  shows  a  plot  pf  the  resulting  error  bound  U$  in  the  azimuthal  angle  estimate 
as  a  function  of  the  true  azimuthal  angle  9,  for  rotation  angle  A 9  =  20°,  when  the  vertical 
arrival  angle  is  <f>  =  20°.  Note  that  the  maximum  error  in  9  occurs  near  the  endfire  directions 
9  =  0°and  9  =  180°.  The  minimum  errors  occur  near  the  broadside  directions  9  =  ±90°. 

5.3.3  Errors  in  Estimated  Vertical  Angle  after  Rotation 

Figure  5.4  shows  the  error  in  the  estimated  vertical  arrival  angle  4>  corresponding  to 
cos  4>±Ueot+  for  the  same  example.  Note  that  the  errors  in  4>  have  the  opposite  trend  to  errors 
in  9.  The  minimum  error  in  vertical  arrival  angle  occurs  near  the  endfire  directions  9  —  0°and 
9  =  180°.  The  maximum  errors  occur  near  the  broadside  directions  9  =  ±90°. 
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Figure  5.2:  Approximate  bound  U$  on  error  in  the  estimated  azimuthal  angle  0  when  the  array 
is  not  rotated,  and  the  vertical  arrival  angle  ^  is  assumed  to  be  zero;  maximum  broadside 
angular  error  is  a0  =  1°;  the  error  h  in  rotation  angle  is  zero;  the  true  value  of  <f>  is  20°. 
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Figure  5.3:  Error  bound  U$  for  the  estimated  azimuthal  angle  6  as  a  function  of  6,  for  vertical 
arrival  angle  <f>  =  20°;  maximum  broadside  angular  error  is  <tq  =  1°;  rotation  angle  A 6  =  20°. 


5.4  General  Trends 


Appendix  A  contains  a  parametric  analysis  of  error  bounds  on  $  and  <f>  resulting  from 
errors  in  the  observations  cos  ,5.  The  parametric  analysis  includes  several  values  of  rotation 
angle  A 9  and  vertical  arrival  angle  <f>. 

The  best  azimuthal  estimates  are  obtained  near  the  broadside  direction;  steep  vertical 
arrival  angles  degrade  the  estimate  of  azimuthal  angle.  If  the  initial  observed  angle  is  near 
broadside  and  the  vertical  arrival  angle  is  very  small,  a  solution  by  rotation  of  the  array  leads  to 
slightly  larger  errors  than  the  simple  uncorrected  estimate.  However,  since  the  vertical  arrival 
angle  is  usually  unknown,  rotation  is  recommended.  Rotation  produces  a  significant  decrease 
in  the  estimated  azimuthal  error,  when  the  vertical  angle  is  steep. 

The  best  vertical  angle  estimates  are  obtained  near  the  endfire  directions;  steep  vertical 
arrival  angles  improve  the  estimate  of  vertical  arrival  angle. 

Appendix  B  contains  a  parametric  analysis  of  r.m.s.  errors  for  the  same  cases.  The  same 
general  trends  apply. 

Appendix  C  contains  a  parametric  analysis  of  the  effect  of  errors  in  the  rotation  an¬ 
gle.  Errors  in  rotation  angle  produce  maximum  errors  in  azimuthal  estimates  near  endfire; 
these  errors  are  independent  of  vertical  arrival  angle.  Maximum  errors  in  vertical  angle  occur 
near  broadside;  these  errors  decrease  as  the  vertical  arrival  angle  increases. 

The  overall  errors  are  minimized  when  the  rotation  angle  is  ±90°.  In  order  to  make 
effective  use  of  the  rotation  solution,  the  rotation  angle  must  be  large  with  respect  to  the 
broadside  resolution  angle  Oq  and  the  rotation  angle  error  h. 


6  Conclusions 


A  Horizontal  Line  Array  (HLA)  is  subject  to  left-right  bearing  ambiguity  due  to  its  one¬ 
dimensional  nature.  It  is  also  subject  to  bearing  bias  error  if  the  signal  vector  is  not  in  the 
same  plane  as  the  HLA.  However,  by  rotating  the  array  to  a  new  orientation,  it  is  possible 
to  estimate  both  the  azimuthal  angle  and  the  vertical  arrival  angle.  The  solution  is  given  in 
closed  form,  for  the  case  of  stationary  source  and  receiver.  The  method  relies  on  the  existence 
of  a  stable  vertical  arrival  angle. 

Graphic  methods  may  be  used  to  aid  in  solution.  However,  a  separate  graph  is  required 
for  each  possible  rotation  angle.  Graphic  methods  are  also  limited  in  accuracy,  since  manual 
interpolation  is  required.  However,  graphs  provide  an  appreciation  of  the  sensitivity  of  the 
solution  to  observation  errors. 

Solution  errors  can  be  expressed  in  closed  form  as  a  function  of  small  errors  in  observed 
angles  and  small  errors  in  the  angle  of  rotation.  Best  estimates  of  vertical  arrival  angle  are  ob¬ 
tained  in  the  end&re  directions,  with  steep  vertical  arrival  angles.  Best  estimates  of  azimuthal 
angle  are  obtained  in  the  broadside  directions,  with  zero  vertical  arrival  angle.  Overall  esti¬ 
mates  improve  as  the  array  rotation  angle  approaches  90°. 

Array  rotation  can  produce  significant  overall  improvement  of  the  azimuthal  angle  esti¬ 
mate.  The  effect  of  large  vertical  arrival  angles  is  reduced  dramatically.  However,  if  the 
vertical  arrival  angle  is  very  small,  angle  estimates  near  broadside  are  slightly  degraded. 
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Appendix  A  Parametric  Analysis  of 
Maximum  Errors 


A  number  of  graphs  have  been  prepared  to  examine  the  sensitivity  of  the  analytic  solution 
to  worst-case  errors  in  the  observed  angles.  The  maximum  positive  and  negative  errors  in  the 
vertical  angle  U and  are  defined  as  follows: 

v;  =  cos“*(C+)  - 

u;  =co!-'(c-)-«  11 


where 


C+  =  max  [(cos^  -  Ucot $),  0] 
C~  =  min  [(cos<*>  +  Ucot^),  1] 


It  is  assumed  in  all  of  these  cases  that  the  rotation  angle  is  known  exactly  (i.e.  h  =  0).  The 
maximum  angle  error  at  broadside  is  ao  =  1°.  Four  rotation  angles  20°,  45°,  90°and  135°are 
considered.  Several  vertical  arrival  angles  are  examined.  Section  5.4  discusses  the  general 
trends  shown  in  these  graphs. 

For  purposes  of  comparison,  figure  A.l  shows  the  maximum  error  Ue  for  the  case  where 
the  array  is  not  rotated,  and  the  vertical  arrival  angle  is  assumed  to  be  zero. 
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Figure  A.l:  Bound  Ue  on  the  maximum  error  in  the  estimated  azimuthal  angle  if  the  array  is 
not  rotated;  a0  =  1°;  h  =  0. 
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Figure  A. 2:  Error  bound  Ut  on  the  estimated  azimuthal  angle  9,  as  a  function  of  6  and  the 
vertical  arrival  angle  <j>,  where  the  maximum  cosine  error  is  sin<r0;  cro  =  1°;  rotation  angle 
A 9  =  20°;  h  =  0. 
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Figure  A. 3:  Error  bounds  Uf  and  U ^  on  the  estimated  vertical  arrival  angle  as  a  function 
of  azimuthal  angle  9,  for  several  vertical  arrival  angles  where  the  maximum  cosine  error  is 
sin<ro;  Co  —  l°i  rotation  angle  A 9  =  20°;  h  =  0. 
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Figure  A.4:  Error  bound  Ug  on  the  estimated  azimuthal  angle  $,  as  a  function  of  0  and  the 
vertical  arrival  angle  <f>,  where  the  maximum  cosine  error  is  sin <70;  oq  =  1°;  rotation  angle 
A9  =  45°;  h  =  0. 
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Figure  A. 5:  Error  bounds  Uf  and  U ^  on  the  estimated  vertical  arrival  angle  4>,  as  a  function 
of  azimuthal  angle  0,  for  several  vertical  arrival  angles  4>,  where  the  maximum  cosine  error  is 
sin<7ol  o'o  =  1°;  rotation  angle  A 9  =  45°;  h  =  0. 
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Figure  A. 6:  Error  bound  Uf  on  the  estimated  azimuthal  angle  0,  as  a  function  of  0  and  the 
vertical  arrival  angle  4>,  where  the  maximum  cosine  error  is  sin oq ;  oq  —  1°;  rotation  angle 
A0  =  90°;  h  =  0. 
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Figure  A. 7:  Error  bounds  Uf  and  U j  on  the  estimated  vertical  arrival  angle  as  a  function 
of  azimuthal  angle  0,  for  several  vertical  arrival  angles  <p,  where  the  maximum  cosine  error  is 
sinaoi  ffo  =  1°;  rotation  angle  A 0  =  90°;  h  =  0. 
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Figure  A. 8:  Error  bound  U$  on  the  estimated  azimuthal  angle  9,  as  a  function  of  9  and  the 
vertical  arrival  angle  where  the  maximum  cosine  error  is  sinao;  oq  =  1°;  rotation  angle 
A 9  =  135°;  h  =  0. 

ROTATION  ANCLE  A 0  -  135.* 
o%  *  1.00*  h  *  0.00* 


AZIMUTHAL  ANGLE  6  IN  DECREES 

Figure  A.9:  Error  bounds  Uf  and  Uj  on  the  estimated  vertical  arrival  angle  as  a  function 
of  azimuthal  angle  9,  for  several  vertical  arrival  angles  where  the  maximum  cosine  error  is 
sin<7o;  <7q  =  1°;  rotation  angle  AS  =  135°;  h  =  0. 


Appendix  B  Parametric  Analysis  of 
R.M.S.  Errors 


A  number  of  graphs  have  been  prepared  to  examine  the  sensitivity  of  the  analytic  solution 
to  random  errors  in  the  observed  angles.  The  positive  and  negative  errors  in  the  vertical  angle 
S+  and  S ^  corresponding  to  the  approximate  r.m.s.  error  Sc0«$  are  defined  as  follows: 

S+ =  cos-1(C+)  -  4> 

S;  =  cos 

where 

C+  =  max  [(cos^>  -  SCOi*),  0] 

C~  =  min  [(cos^  +  Sco,$),  1] 

It  is  assumed  in  all  of  these  cases  that  the  rotation  angle  is  known  exactly  (i.e.  h  = 
0).  The  broadside  resolution  angle  oq  is  1°.  Four  rotation  angles  20°,  45°,  90°and  135°are 
considered.  Several  vertical  arrival  angles  are  examined.  Section  5.4  discusses  the  general 
trends  shown  in  these  graphs. 

For  purposes  of  comparison,  figure  B.I  Bhows  the  r.m.s.  error  Se  for  the  case  where  the 
array  is  not  rotated,  and  the  vertical  arrival  angle  is  assumed  to  be  zero.  For  this  special  case, 
the  approximate  r.m.s.  errors  are  defined  as  follows: 


(B  2) 


(B.I) 


S|«  [BJ(M)  +  A£2(0,<A,O  +  >>2]*  ;  A0  =  o 
S*  =  <(>  ;  AO  —  0  . 
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Figure  B.l:  Approximate  r.m.s.  error  Sg  for  the  estimated  azimuthal  angle  if  the  array  is  not 
rotated;  oq  =  1°;  h  =  0. 
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Figure  B.2:  Approximate  r.m.s.  error  S#  in  azimuthal  angle  9,  as  a  function  of  9  and  the  vertical 
arrival  angle  <f>,  where  the  r.m.s.  cosine  error  is  sin  <tq;  ffo  =  1°;  rotation  angle  A 9  =  20°;  h  =  0. 
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Figure  B.3:  Errors  Sf  and  S j  in  vertical  arrival  angle  as  a  function  of  azimuthal  angle  9, 
for  several  vertical  arrival  angles  where  the  r.m.s.  cosine  error  is  sin  ct0;  &o  =  1°;  rotation 
angle  A 9  =  20°;  h  =  0. 
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Figure  B.4:  Approximate  r.m.s.  error  S$  in  azimuthal  angle  as  a  function  of  6  and  the  vertical 
arrival  angle  <f>,  where  the  r.m.s.  cosine  error  is  sin  er0;  <r0  =  1°;  rotation  angle  A 9  =  45°;  h  =  0. 
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Figure  B.5:  Errors  S£  and  S ^  in  vertical  arrival  angle  <f>,  as  a  function  of  azimuthal  angle  8, 
for  several  vertical  arrival  angles  where  the  r.m.s.  cosine  error  is  sinaoi  &o  —  1°;  rotation 
angle  A 9  =  45°;  h  =  0. 
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Figure  B.6:  Approximate  r.m.s.  error  Sg  in  azimuthal  angle  0,  as  a  function  of  0  and  the  vertical 
arrival  angle  <)>,  where  the  r.m.s.  cosine  error  is  sin  <ro;  er0  =  1°;  rotation  angle  A0  =  90°;  h  =  0. 
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Figure  B.7:  Errors  Sf  and  Sj  in  vertical  arrival  angle  <f>,  as  a  function  of  azimuthal  angle  6, 
for  several  vertical  arrival  angles  4>,  where  the  r.m.s.  cosine  error  is  sincro;  ctq  =  1°;  rotation 
angle  AO  =  90°;  h  =  0. 
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Figure  B.8:  Approximate  r.m.s.  error  S«  in  azimuthal  angle  9,  as  a  function  of  9  and  the  vertical 
arrival  angle  <f>,  where  the  r.m.s.  cosine  error  is  sin^ol  =  1°;  rotation  angle  A 9  =  135°; 


h  =  0. 
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Figure  B.9:  Errors  Sf  and  S ^  in  vertical  arrival  angle  <f>,  as  a  function  of  azimuthal  angle  9, 
for  several  vertical  arrival  angles  <f>,  where  the  r.m.s.  cosine  error  is  sinco;  <^o  =  1°;  rotation 
angl^  A 9  =  135°;  h  =  0. 


Appendix  C  The  Effect  of  Error  in 
Rotation  Angle 


A  number  of  graphs  have  been  prepared  to  examine  the  sensitivity  of  the  analytic  solution 
to  an  error  h  in  the  rotation  angle  A 0.  In  all  of  the  examples,  h  =  1°. 

It  is  assumed  in  all  of  these  cases  that  the  observations  cos^i  and  cos  Pi  are  exact,  i.e. 
cr0  =  0°.  Four  rotation  angles  20°,  45°,  90°and  135°are  considered.  Several  vertical  arrival 
angles  are  examined.  Figures  C.2  to  C.9  show  the  maximum  errors  resulting  from  a  maximum 
error  of  l°in  the  rotation  angle.  Figures  C.ll  to  C.18  show  the  r.m.s.  errors  resulting  from  an 
r.m.s.  error  of  1°. 

For  purposes  of  comparison,  figures  C.l  and  C.10  show  the  maximum  error  Ug  and  the 
r.m.s.  error  Sg  for  the  case  where  the  array  is  not  rotated,  and  the  vertical  arrival  angle  is 
assumed  to  be  zero. 


C.l  Maximum  Errors  Due  to  Error  in  Rotation  Angle 
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F:c  :re  C.l:  Bound  £/#  on  the  maximum  error  in  the  estimated  azimuthal  angle  if  the  array 

not  rotated;  <tq  =  0°;  h  =  1°. 
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Figure  C.2:  Error  bound  Ug  on  the  estimated  azimuthal  angle  0,  as  a  function  of  0  and  the 
vertical  arrival  angle  <t>,  where  the  maximum  cosine  error  is  Binooi  <ro  =  0°;  rotation  angle 
A0  =  20  °;h=  1°. 
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Figure  C.3:  Error  bounds  Ut  and  on  the  estimated  vertical  arrival  angle  as  a  function 
of  azimuthal  angle  0,  for  several  vertical  arrival  angles  where  the  maximum  cosine  error  is 
sin^o;  <tq  =  0°;  rotation  angle  A0  =  20°;  h  =  1°. 
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Figure  C.4:  Error  bound  Ug  on  the  estimated  azimuthal  angle  0,  as  a  function  of  6  and  the 
vertical  arrival  angle  4>,  where  the  maximum  cosine  error  is  sinao;  <^o  =  0°;  rotation  angle 
AO  =  45°;  h  =  1°. 
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Figure  C.5:  Error  bounds  Uf  and  U j  on  the  estimated  vertical  arrival  angle  as  a  function 
of  azimuthal  angle  0,  for  several  vertical  arrival  angles  4>,  where  the  maximum  cosine  error  is 
sincro;  ffo  =  0°;  rotation  angle  A0  =  45°;  h  =  1°. 
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Figure  C.6:  Error  bound  Ug  on  the  estimated  azimuthal  angle  0,  as  a  function  of  6  and  the 
vertical  arrival  angle  <f>,  where  the  maximum  cosine  error  »  sin<70;  <to  =  0°;  rotation  angle 
A0  =  90°;  h  =  1°. 
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Figure  C.7:  Error  bounds  Uf  and  Uj  on  the  estimated  vertical  arrival  angle  4>,  as  a  function 
of  azimuthal  angle  0,  for  several  vertical  arrival  angles  where  the  maximum  cosine  error  is 
sin<70;  0o  =  0°;  rotation  angle  A0  =  90°;  h  =  1°. 
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Figure  C.8:  Error  bound  Ug  on  the  estimated  azimuthal  angle  0,  as  a  function  of  0  and  the 
vertical  arrival  angle  <f>,  where  the  maximum  cosine  error  is  Binaol  <ro  =  0°;  rotation  angle 
A0  =  135°;  h  =  1°. 
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Figure  C.9:  Error  bounds  Uf  and  U ^  on  the  estimated  vertical  arrival  angle  <t>,  as  a  function 
of  azimuthal  angle  6,  for  several  vertical  arrival  angles  where  the  maximum  cosine  error  is 
sinao;  =  0°;  rotation  angle  A0  =  135°;  h  =  1°. 
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C.2  R.M.S.  Errors  Due  to  Error  in  Rotation  Angle 
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Figure  C.10:  Approximate  r.m.s.  error  Sg  for  the  estimated  azimuthal  angle  if  the  array  is  not 
rotated;  <7q  =  0°;  h  —  1°. 


43 


ROTATION  ANCLE  LB  -  20.* 
o,  =  0.00*  h  -  1.00* 


VERTICAL  ANCLE  4> 

-  0* 

-  -  •  5.* 

- 10.’ 

—  -20.* 

- SO.* 

40.* 

- 50.’ 


Figure  C.ll:  Approximate  r.m.s.  error  S$  in  azimuthal  angle  8,  as  a  function  of  8  and  the 
vertical  arrival  angle  <?,  where  the  r.m.s.  cosine  error  is  sinoo;  oq  —  0°;  rotation  angle 
A 8  =  20°;  h  =  1°. 
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Figure  C.12.  Errors  Sf  and  S ^  in  vertical  arrival  angle  as  a  function  of  azimuthal  angle  8 , 
for  several  vertical  arrival  angles  <f>,  where  the  r.m.s.  cosine  error  is  sin^o;  <*o  =  0°;  rotation 
angle  A 9  —  20°;  h  =  1°, 
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Figure  C.13:  Approximate  r.m.s.  error  Se  in  azimuthal  angle  0,  as  a  function  of  6  and  the 
vertical  arrival  angle  <j>,  where  the  r.m.s.  cosine  error  is  sinoo',  <r0  =  0°;  rotation  angle 
AN  =  45°;  h  =  1°. 
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Figure  C.14:  Errors  Sf  and  S ^  in  vertical  arrival  angle  as  a  function  of  azimuthal  angle  $, 
for  several  vertical  arrival  angles  4>,  where  the  r.m.s.  cosine  error  is  sin  (To;  <*o  =  0°;  rotation 
angle  AN  =  45°;  h  =  1°. 
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Figure  C.15:  Approximate  r.m.s.  error  St  in  azimuthal  angle  6,  as  a  function  of  6  and  the 
vertical  arrival  angle  <j>,  where  the  r.m.s.  cosine  error  is  sinooi  &o  =  0°;  rotation  angle 
A6  =  90°;  h  =  1°. 
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Figure  C.16:  Errors  Sf  and  S ^  in  vertical  arrival  angle  as  a  function  of  azimuthal  angle  9, 
for  several  vertical  arrival  angles  where  the  r.m.B.  cosine  error  is  sin<r0;  ffo  =  0°;  rotation 
angle  A8  =  90°;  h  =  !°. 
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Figure  C.17:  Approximate  r.m.s.  error  Sg  in  azimuthal  angle  9,  as  a  function  of  9  and  the 
vertical  arrival  angle  <f>,  where  the  r.m.s.  cosine  error  is  sino0;  ffo  =  0°;  rotation  angle 
A9  =  135°;  h  =  1°. 
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Figure  C.18:  Errors  Sf  and  in  vertical  arrival  angle  <f>,  as  a  function  of  azimuthal  angle  9, 
for  several  vertical  arrival  angles  4>,  where  the  r.m.s.  cosine  error  is  sincr0;  oq  =  0°;  rotation 
angle  A 9  =  135°;  h  =  1°. 
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be  expressed  in  closed  form. 
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